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ABSTRACT 

We examine atmospheric heating by radio active galactic nuclei (AGN) in distant X-ray clusters by cross 
correlating clusters selected from the 400 Square Degree (400SD) X-ray Cluster survey with radio sources 
in the NRAO VLA Sky Survey. Roughly 30% of the clusters show radio emission above a flux threshold of 
3 mJy within a projected radius of 250 kpc. The radio emission is presumably associated with the brightest 
cluster galaxy. The mechanical jet power for each radio source was determined using scaling relations between 
radio power and cavity (mechanical) power determined for nearby clusters, groups, and galaxies with hot 
atmospheres containing X-ray cavities. The average jet power of the central radio AGN is approximately 
2 X lO'^^ergs"^. We find no significant correlation between radio power, hence mechanical jet power, and 
the X-ray luminosities of clusters in the redshift range 0.1 - 0.6. This implies that the mechanical heating rate 
per particle is higher in lower mass, lower X-ray luminosity clusters. The jet power averaged over the sample 
corresponds to an atmospheric heating of approximately 0.2 keV per particle within Rsoo- Assuming the current 
AGN heating rate does not evolve but remains constant to redshifts of 2, the heating rate per particle would 
rise by a factor of two. We find that the energy injected from radio AGN contribute substantially to the excess 
entropy in hot atmospheres needed to break self-similarity in cluster scaling relations. The detection frequency 
of radio AGN is inconsistent with the presence of strong cooling flows in 400SD clusters, but does not exclude 
weak cooling flows. It is unclear whether central AGN in 400SD clusters are maintained by feedback at the 
base of a cooling flow. Atmospheric heating by radio AGN may retard the development of strong cooling flows 
at early epochs. 

Subject headings: Galaxies: clusters: general; Galaxies: clusters: intracluster medium; Galaxies: quasars: 
general; X-rays: galaxies: clusters; Radio continuum: galaxies 
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1. INTRODUCTION 

Heating by active galactic nuclei (AGN) may be re- 
sponsible for two significant phenomena in galaxy clus- 
ters. First, cool, dense cores found in the X-ray atmo- 
spheres of clusters are expected to cool efficiently to low 
temperatures (F abian 199 4). However, the levels of cold gas 
and star formation lie well below the values expected from 
pure radiative cooling (reviewed by lPeterson & Fabianll2006t 
iMcNamara & Nulsen 2007). The generally accepted solu- 
tion is that the cooling is being compensated by one or more 
heating mechanisms, possibly including conduction from the 
hot g a s in non-cooling regions (e.g . jNaravan & Medvedev 
200U IZakamska & NaravSil 12003^ IVoigt & Fabian 2007 



Wise et alJ 12004. and AGN h eating (e.g. iB innev & Tab^ 
1995'; 'Churazov et all l2002t iRoychowdhury et al.i i2004 
Voit & Donahue 20051)- Among the possible heating mecha- 
nisms, AGN heating stands out because of both the enormous 
energy released and the opportunity for feedback, since AGN 
power is controlled by the accretion rate. Second, galaxy 
groups and poor clusters are less luminous than expected 
based on the Lx -T relation derived for massive clusters 
jMarkevitch" 1 998"; 'Arnaud & Evrard"1999'; 'Ruszk owski et all 
|2Q04; San derson et al. 2003 ; Popesso et a l. 2005). This is in- 
terpreted as a break of the self-similar scaling relation be- 
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tween X-ray luminosity and gas temperature expected in 
structure form ation models that i nclude only gravitational 
heating (iKaiserill986tlEvrard et"al][l996.) . The lower than ex- 
pected X-ray luminosity for a given temperature is thought 
to reflect extra entropy (heat) in the gas that was injected 
by star formation and AGN at early tim es (e.g. Kaiser 1991; 
Balogh et al. 1999"; 'Borgani et al.' 2005; Ponman et al. 2003; 
Croston et al. 2005; Sanderson et al. 2005; Jetha et al. 2007 



Cavagnolo et al.l l2()09l) ! The "preheating" model proposes 



that the excess entropy wa s injected during the early stages 
of structure formation (e.g. lKaiseiJll991l; iBalogh et al.lll999l; 
iMcCarthv et al.ll200l . 

The amount of energy released by an AGN, on the or- 
der of 10^^ AfBH/[10^ Mq] erg, is more than enough to be 
the major heating source. Nevertheless, how and where 
the energy is d istributed into the ICM rem ains an issue. 
As discussed by IMcNamara & NulsenI (I2007h . AGN heat in 
two ways. First, AGN outbursts inflate cavities. The en- 
thalpy released as the cavities rise buoyan tly is dispersed 
into the ICM (B ruggen & Kaiser 2002; Revnolds et al.l l2002l; 
IChurazov et al.l 12002; Birzan et al. 2004^! Therm alization 
of the energy of cosmic rays diffusing out of a radio lobe 
is an alter native channel for converting lobe free energy to 
heat (e.g. [Mathews '2009). Second, the AGN inject en- 
ergy into the IC M through shocks (e.g . Forman et al. 20051 
IMcNamara et al. 2005; N ulsen et al.ll200 7) and sound waves 
(e.g.l Fabian et al...200 5a). Cavity power can be estimated us- 
ing the enthalpy contained in the cavity and its rise time. In 
practice, the enthalpy is calculated as 4pV^av, assuming the 
cavity is filled with relativistic particles, and the age can be 
estimated using the terminal velocity of the buoyantly rising 
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bubbles (iBirza n et al.ll2004l: iDunnet aljllool iRaffertv et ai] 
[2006; Birzan et al.ll2008h . All of these quantities, the pressure 
of the surrounding gas p, the volume of the cavity Vcav and 
sound speed in the ICM, c an be obtained from X- ray data in 
principle. Simulations by (Mendy gral et aljr201 Ih have veri- 
fied that the cavity power estimated using this approach can 
be accurate to approxim ately a factor of two. Furthermore, 
IBirzan et al.1 (12008,) and .Cavagnoloet al.l (120101) find that the 
radio power and the buoyant cavity power are correlated for 
samples ranging from elliptical galaxies and galaxy groups to 
massive galaxy clusters. Thus, the radio power can be used 
as a proxy for the cavity power for large samples of radio 
galaxies embedded in hot atmospheres that lack deep X-ray 
observations. However, cavity power estimated using this ap- 
proach is probably a lower limit, because it does not include 
the additional energy released by shocks (Forman et al. 2005; 
iMcNamara et al .''20Q5: Nulsen et al. 2005a b). 

It is more challenging to calculate the total AGN energy 
deposited over the life of a cluster Since the power output of 
AGN is highly variable and outbursts are intermittent, current 
AGN power is not necessarily representative of its historical 
average. Even under the simplest assumption, that the power 
of an AGN is the same in eve ry outburs t (which is incon sis- 
tent with data, e.g. lFabian eFaJ . 2006; R andall etal.ll20 111) , an 
estimate is still required for the fraction of the time that the 
AGN is in outburst, i.e., for its duty cycle, which is also inac- 
cessible for individual AGN. An alternative to studying single 
systems is to determine the average AGN power for a large 
sample. For a well chosen sample, the average power should 
represent the time average of the AGN power for individual 
sample members. The key to success here is to use a well- 
defined, unbiased, and sufficiently large sample. 

In this paper, we estimate the average jet power in distant. 
X-ray selected clusters found in the 400 Square Degree Clus- 
ter Survey (400SD; Burenin et al. 2007; Vikhlinin et al. 2005|) 
by cross coiTelating t he 400SD catalog w ith the NRAO VLA 
Sky Survey (NVSS; ICondon et all 11998.) . We calculate jet 
powers for radio sources above a flux hmit of '^ 3 mJy. Fol- 
lowing earlier work on AGN heating in cooling core clusters 
dPunn & Fabian 2006, 2008; Rafferty et al. 2006, 2008), we 
evaluate the AGN power injected into normal clusters to ad- 
dress the question of how much energy is deposited by AGN 
into cluster halos over a significant fraction of their ages. Note 
that AGN feedback would be expected to enhance the inci- 
dence of radio outbursts at the centers of cooling core clusters. 
As a result, this question cannot be answered based on a sam- 
ple of nearby cooling core clusters, which does not represent 
the general cluster population. Also, since we are interested in 
the total power injected by AGN, multiple radio sources may 
be attributed to a single cluster in the present study. 

The paper is arranged as follows. In Sj2] we introduce the 
400SD clusters, and the matched radio sources. In pa rticular, 
we discuss the background source correction in Sj23] and ex- 
amine the morphology of a subsample of the NVSS sources 
using images from the su rvey Faint Imag es of the R adio Sky 
at Twenty-cm (FIRST; iBecker et al.l[T99 5) in ^g^ Then we 
present our analysis and discuss the interpretation: the corre- 
lation between the 1 .4 GHz power of the radio AGNs and the 
X-ray luminosity of their host clusters in i l3.1l the radial dis- 
tribution of radio sources in the clusters in i l3.2l the jet power 
estimated from the radio power in 0.3I and the estimation of 
AGN heating energy in i i3.4l We summarize our results in ^ 



2. SAMPLE 



2.1. X-ray data 

The 400SD survey is one of the largest serendipitous X-ray 
surveys of clusters based on ROSAT Position Sensitive Pro- 
portional Counter pointed observations. It includes 242 op- 
tically confirmed clusters and groups brighter than an X-ray 
flux Hmit of 1.4 x 10^^-^ ergs^^cm^^ in the 0.5 — 2keV en- 
ergy band, and lying at redshifts z < 0.9 within a 397 deg^ 
area. It is arguably one of the best surveys for studying the 
evolution of clusters because the 400SD samples "normal" 
clusters at higher redshift. It is not targeted at the most mas- 
sive clusters. It includes clusters with Lx,boi ~ 10'*'* ergs~* 
at z > 0.3 (Fig. [D, in contras t to the all-sky surveys, such 
as MACS (Ebelinge t all 1200 Ih . BCS-neBCS ("E behng et all 
1997, 2000), which are targeted at the rarer, more luminous 
and more massive clusters. In addition, the ratio of strong 
cooling core clusters (SCC ; defined as cooling time less than 
1 Gvr in iMittal et al J 120091) in 400SD is very low according 
to ISantos et all (120101) . who find none SCC out of 20 clus- 
ters at z > 0.5), comparing to the all-sky surveys, e.g. BCS 
(|DuniL& Fabian'2008) and HIFLUGCS (iMittal et aklifOOl 
A similar conclusion is also made by S amuele et alJ (120111) 
studying emission lines in the optical spectra of the brightest 
cluster galaxies (BCGs) in 400SD clusters. 

Here we provide a short summar y of the measuremen ts of 
the 400SD clusters in the catalog of lBurenin et alj (12007!) and 
Vikhlinin etal. (1998). 

• The clusters were confirmed optically using images 
with magnitude limit mj^i ~ 24. Their redshifts were 
measured spectroscopically. 

• The center position and core radius of each cluster were 
derived from the best-fit /3-model with (3 fixed to 0.67. 
The typical systematic error of the cluster positions is 
^ 17" with respect to the brightest cluster galaxy. 

• The X-ray flux in the energy band 0.5 — 2.0 keV was 
estimated as the mean of two fluxes calculated from 
the integration of best-fit /3-models with j3 fixed to 0.6 
and 0.7 respectively. The 10% systematic uncertainty 
caused by the uncertainty of /3 value was added. 

• The X-ray Luminosity Lx, 0.5-2.0 was derived interac- 
tive ly assuming the Lx-T r elation ofM arkev itch (19981) 
and lFukazawa et al.1 (119981) . 

• To be consistent, we extrapolated to obtain the X- 
ray bolometric luminosity Lx.boi (Fig-ID assuming the 
same Lx-T relation. 

We used the high resolution ACIS/Chandra data, which 
is available for a subsample of the 400SD clusters 
(Vikhlinin et al. 2009, 52 clusters in total), to examine the un- 
certainties in the luminosities and centroids determined from 
the PSPC data. Compared to centers determined visually from 
the ACIS images (see Fig.|2|, we find that the average offset 
of the position in the catalog is ^ 12", which is consistent 
with the uncertainty quoted above. In the rest frames of these 
clusters, all but one of the offsets correspond to projected dis- 
tances of less than ^ 250kpc. The exception is a binary sys- 
tem (CLJ0152. 7-1357) whose cataloged centroid is located 
bet ween the two clusters which are separated by ~ 100". 

In lBurenin et alj (120071) . the uncertainty in the luminosities 
derived from the PSPC fluxes using more accurate Chandra 
fluxes is 19 ± 0.6%. In addition, a few outliers, which are 
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Fig. 1 . — X-ray bolometric luminosity distribution with redshift. Red dots 
mark the clusters with matched radio sources within 2 Mpc radius. Those 
clusters with radio sources at the center (r< 250 kpc) are plotted with blue 
circles. For the clusters located at the region of NVSS survey and no radio 
sources detected are marked with black spheres. The other 400SD clusters 
with z < 0.1 are marked with dots. The Y-axis is Bolometric luminosity 
derived using PSPC flux in 0.5 — 2.0 keV assuming Lx-T relation (see JTTJ. 



likely caused by conta mination from X-ray point sources, can 
be found in Fig. 23 of iBurenin efal] ( 120071) . 



in projection and are not physically associated with the clus- 
ters. We discuss the background correction in ; 
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Fig. 3 . — Radio luminosity of NVSS sources around clusters. The luminos- 
ity density of matched NVSS sources are plotted with respect to the redshift 
of their associated cluster. The red spheres mark the sources within 250 kpc 
radius of cluster core, and the smaller black dots mark those within 2 Mpc 
radius. 
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Fig. 2. — Offsets between cluster centers measured using ACTS and PSPC 
data. The offsets are estimated from a subset (52) of clusters with ACTS im- 
ages. The upper histogram shows angular separation and the lower histogram 
shows projected distance at the cluster redshift. The mean offset is marked 
by a red vertical line. 



2.2. NVSS sources 

In order to match the coordi nates of 400Sp clus ters with 
1.4Ghz sources in the NVSS JCondonet al.lll99 8'). we ex- 
tracted a subsample consisting of the 400SD clusters in the 
region covered by the NVSS survey. This procedure yields 
196 clusters at z > 0.1. We find, in total, 782 radio sources 
above a flux limit of 3 mJy 0, located within a projected dis- 
tance of 2 Mpc of 166 clusters at z > 0.1. Of these sources, 61 
are located within a 250 kpc radius of 56 clusters at z > 0.1. 
Consideri ng t he uncertainty in the centroids for the 400SD 
clusters (i )2.1l i and the resolution of NVSS images, all sources 
within 250 kpc are consistent with being associated with the 
central galaxies of the clusters. We note that a small fraction 
of these matched radio sources are background sources seen 



^ The completeness reaches 90% at 3 mJy jCondon et alj|199"8l) . We also 
tried a lower flux limit, 2 mJv, and found that the results are qualitatively 



insensitive to the limit (see ^3.4 
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Fig. 4. — Radial distribution of the surface density of NVSS sources 
around 400SD clusters. The NVSS sources are split into two redshift ranges: 
0.1 < z < 0.3 (black) and z > 0.3 (red dashed) histograms. The gray 
histogram shows the surface density distribution for the low redshift sample 
with a luminosity cut P1.4 > 1 X 10^"' WHz^^. 

We calculated the radio power density. Pi 4 (Fig. O, for 
each matched NVSS source, using the cluster redshift and a 
spectral ind ex of —0.8, which is the average for NVSS sources 
JCondon et al. 1998). The cluster redshift was used because 
we cannot identify and measure the redshift of the host galaxy. 

In Fig. m we plot the surface density of detected NVSS 
sources around the 400SD clusters as functions of distance 
from the X-ray centroid, for the redshift bins, z > 0.3 and 
0.1 < z < 0.3. The surface density of NVSS sources as- 
sociated with clusters declines more quickly with distance at 
higher redshifts (red) compare to those at redshift between 0. 1 
and 0.3 (black). The distribution for the lower redshift sample 
is flat beyond ^ 250 kpc. The difference between the profiles 
arises because of the radio flux limit. Our adopted flux limit of 
3 mJy corresponds to a minimum luminosity of lO^"* W Hz^^ 
(see Fig. [3]) for radio sources beyond z — 0.3, and drops to 
10^3 WHz~^ for radio sources at the redshift 0.1. Thus, we 
detect many more low power radio sources within and pro- 
jected onto the nearby clusters compared to those at higher 
redshifts. Plotting the radial distribution of only the high lu- 
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minosity sources (Pi. 4 > 10^"* WHz^^) for clusters lying be- 
tween 0.1 < z < 0.3 (gray histogram in Fig.|4|i gives a closer 
match to the distribution for clusters beyond z = 0.3. How- 
ever, the surface density of NVSS sources within the central 
250 kpc is significantly higher in the higher redshift clusters 
(red dashed histogram) compared to the lower redshift clus- 
ters (gray histogram). 

2.3. Background correction 
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Fig. 5. — Probability that a detected source is from the background. Pro- 
jected distances for the sources located within 250 kpc are set to 250 kpc, 
reflecting the uncertainty in the cluster centers. 
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Fig. 6. — FIRST images of matched NVSS sources. We show the six 
strongest sources found within 250 kpc, and the images ( 3' X 3' ) are sorted by 
decreasingflux, from top left to bottom right. The coordinates of the sources 
are noted on the top of each panel. 



We are unable to identify the radio host galaxies due to the 
poor resolution of the NVSS survey. Furthermore, we do not 
have their redshifts. We have therefore corrected our numbers 
for background contamination statistically. We assumed that 
the spatial distribution is locally homogeneous and the flux 
function is independent to the location on the sky. Thus, the 
probability of finding any background sources within an area 
A can be expressed as 



Pbkg(A) = 1 



-A 



where 



A = 



^bkg 



A 



bkg 



xA (1) 



the mean A is estimated from the number density of the back- 
ground sources, the number of the NVSS sources in an an- 
nulus from 2° to 4° around the cluster of interest is A^bkg, 



and ^bkg is the area of this annulufl The typical number 
of sources in the annulus is about 1600 and the background 
number density is about 42 deg~^, consistent with the num- 
ber density mentioned in ICondon et al. (1998). Fig.|5]shows 
that any detected source is likely to be a background source 
at large radii. Therefore, we confine out attention mostly to 
sources located within 250 kpc of the center 

In order to correct the flux contamination, we calculated the 
expected flux of background sources within the background 
area defined above. When we calculated the average flux of 
radio sources within some projected distance of a cluster, we 
scaled the expected background flux according to the ratio of 
area, and subtracted it from the gross average flux. 

2.4. High resolution radio source images for a subset of the 

sample 

Before we investigate the statistical properties of the match- 
ing radio sources, we briefly discuss the properties of a few 
sources with higher resolution radio images obtained from the 
FIRST data archive. Overall, more than half of the NVSS 
sources were observed in the FIRST survey. For example, we 
examined the images of 33 out of 56 NVSS sources matched 
within a 250 kpc radius of the 400SD clusters. We find 14 of 
these sources have double radio source morphologies, and 7 
of the remaining 19 single sources are resolved at the FIRST 
resolution (^ 5") with a minimum FIRST flux ^ 1.9 mJy . We 
also checked the optical images of Vikhlinin et alj (12009 ), and 
find that ^ 75% (25/33) of the radio sources within 250 kpc 
are associated with BCGs. We are unable to identify the 
BCGs of the remaining 8 clusters with radio sources using 
the available optical imagery. 

In Fig. |6] we show FIRST images for the strongest six 
sources located within 250 kpc radius of the 400SD clusters. 
Of particular interest is 3C288, which is the brightest NVSS 
source in our sample (Pi 4 = 700 x lO^^'WHz"^). Its 
host is the BCG of cluster 400DJ1338H-3851 at z = 0.246 
with Lx, 0.5-2.0 = 9.6 X 10^'^ ergs^^. It is a well-studied 
object b ecause of its bright but edge-darkened radio mor- 
phology ([Bridle et alJ[T98^ISaripam & Subrahmanvanl2009l: 



iLal et al.ll2010l) , and ha s been characterized as a transitional 
FRI/FRII radio galaxy (Fanar off & RilevllI974l) . In the core 
of the cluster, Lai et al. (2010) detected two surface bright- 
ness discontinuities, which are interpreted as shocks induced 
by the supersonic inflation of the radio lobes. The inferred 
shock energy is comparable to the cavity enthalpy of Pjot = 
6.1 X 10^'' ergs-i. 



3. RESULTS AND DISCUSSION 



3.1. 



Correlation between radio power ofAGN and X-ray 
luminosity of clusters? 

Using our catalog of 400SD clusters, we now examine the 
relationship between the radio power of the central galax- 
ies and cluster X-ray luminosity. A simple flux-flux plot 
in Fig. I2] shows no correlation between the radio fluxefl 
(Fi 4) of the sources within 250 kpc radius and the X-ray 
fluxes (Fx,o.5-2.okcv)- The Spearman correlation is only 
0.15 ±0.20. To make sure the correlation is not hidden by red- 
shift dependent factors, such as reddening, we plot the ratio of 

* We have also tried the background source density calculated using the 
entire NVSS source catalog, and found no significant difference in our results. 

' Radio fluxes are summed for all sources within 250 kpc of each cluster 
cener However, most of the clusters (51 out of 56) with radio sources within 
250 kpc radius only contain a single radio source, so that the fluxes could be 
taken as those of individual radio sources. 
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radio power to X-ray luminosity versus redshift in Fig.|8] Al- 
though the points seem to be centered on a constant ratio of 
^ 0.001, the lower limits, estimated for clusters with no de- 
tected NVSS sources (Fi 4 > 3 mJy), and the upper limit (red 
curve) for clusters not detected at the X-ray flux limit of the 
400SD survey, suggest that any correlation is deceptive. In 
summary, no correlation is found between radio power and 

X-ray luminosity in our cluster sample. 

Several previous studies (e.g. Sun 200"9t iMittal et al.ll2009h 
have demonstrated that the correlation between X-ray lumi- 
nosity of the cooling core and the radio power only appears 
in "strong" cooling core clusters. Therefore, the lack of a 
correlation between radio and X-ray power is consistent with 
the l ack of strong cooling core clusters among 400SD clus- 
ters dVikhUnin et al.ll200^ ISantos et"ani2010t ISamuele et all 



120111) . Of the52 400SD clusters with Chandra data, 20 have 
NVSS sources within 250 kpc of the center. Of those, 6 
clusters have surfa ce brightness concentrations csb > 0.1 
(ISantos et al.l 120081) . c orresponding to c entral cooling times 
shorter than ~ 5 Gy (ISantos etal.1 1201 Oi). These are likely 
to be weak cool core clusters (WCC in iMittal et al] 120091) . 
These sources are plotted in red in Fig. Q and their radio and 
X-ray fluxes do appear better correlated than for the bulk of 
the sample, consistent with previous studies. 

We also consider the dependence of the fraction of clusters 
with NVSS radio sources on X-ray luminosity. In Fig.|9] the 
incidence of radio sources appears to increase slightly with 
X-ray luminosity. However, the range of X-ray luminosi- 
ties sampled depends on the redshift in a flux-limited survey. 
The low luminosity clusters in the sample are only detectable 
at lower redshift because of the X-ray flux limit. Separat- 
ing the sample into two redshift slices (0.1 < z < 0.3 and 
0.3 < z < 0.5 in Fig. |9]l, we find that the apparent corre- 
lation between X-ray luminosity and the fraction of cluster 
with NVSS sources disappears. Thus, we suggest that the dif- 
ference between the blue points (higher redshift, more lumi- 
nous clusters) and the red points (lower redshift, less luminous 
clusters ) is probab l y caus ed by redshift and not X-ray lumi- 
nosity. iBest et ai] ( 120071) found no correlation between the 
fraction of radio-loud AGN in cluster galaxies and cluster ve- 
locity dispersion, which is, a proxy f or X-r ay luminosity and 
cluster mass. However, L in & Mohj (l200 7h found a clear in- 
crease in the "radio active fraction" in a massive cluster sam- 
ple selected when they have bright BCGs with M^ < —24. 
In order to examine this apparent inconsistency, we need to 
study the k-band luminosity of the BCGs which is beyond the 
scope of this paper 

Although we see no evidence of a dependence on X- 
ray luminosity, the incidence of radio AGN does in- 
crease marginally at higher redshift. This may reflect 
higher AGN fractions in the past, as suggested by some, 
ej.,by Martini etal. (2007, 2009) for X-ray AGN and by 
iGalametz et al.l (2009) for AGN selected by radio, infrared 
and X-ray. However. Gralla et al..(.2010.) analyzed radio AGN 
in the 618 clusters of the Red Sequence Clusters Survey at 
redshifts 0.35 < z < 0.95, finding that the average number of 
AGN is independent of redshift. Although they are less sig- 
nificant, our results agree better with the trend found for radio 
AGN by IGalametz et al.l (2009), of an increase in incidence 
by about a factor of 2 (~ 2a significance) from z < 0.5 to 
0.5<z< 1.0. 

3.2. Radial Distribution and Incidence of radio AGN in 
clusters 



In Fig. [TOl we plot the average cumulative number of ra- 
dio AGN within radius r versus r for the 400SD clusters. 
The number of radio sources increases from an average of 
^ 0.3 in the cluster cores to ^ 1 within 2 Mpc. However, 
the number of radio sources increases much more slowly than 
the area (i.e., ex r^), consistent with the observation that 
the surface density of radio sources peaks at cluster centers 
(iLedlow & Owenlll99l [Best et alJl2007t iLin & MohriilOOTJ: 
lCroftetal.ll2007iK 

Since most 400SD clusters have no more than one radio 
source within 250 kpc, the average number of radio sources 
shown as the leftmost point in Fig.[TO]is equivalent to the ra- 
dio source fractions in Fig. |9l The average number of radio 
sources we find in the cores of 400SD clusters is significantly 
lower than the fraction of radio-loud BCGs in nearby cooling 
core clusters, which li es between 75% to 100% (Burns 199(1 
lEdwards et all 120071: iDunn et al.l HOO^; Mit tal et alJ 12009). 
This fraction is expected to correlate with the central cool- 
ing time. Our detection fraction is close to the fraction of ra- 
dio sources f ound in nearby, non -cooling core clusters (45%) 
as defined bv IMittal et alj (120091) . and is identical to the frac- 
tion found in optically selected clusters (30%) determined by 
Best et al. (2007). While a more quantitative comparison be- 
tween the numbers is difficult, a low radio source fraction for 
the 400SD clusters is consistent with our claim that most of 
them lack strong cooling cores. At the same time, our detec- 
tion fraction is significantly higher than the radio loud fraction 
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Fig. 7. — F1.4 versus Fx for the 400SD clu sters with NVSS sources within 
250 kpc. The expected background fluxes ([(23) have been subtracted, al- 
though the correction is insignificant for this small region. Using the archived 
Chandra data, we can estimate the cooling time of clusters referring from the 
concentration parameter csb- The clusters with cooling time less than 5 Gyr 
are plotted in red. 
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Fig. 8. — Ratio of radio power within 250 kpc to X-ray luminosity for 
400SD clusters. Dots mark clusters with detected radio sources and upper 
hmits are for non-detections, assuming one 3 mJy radio source in each cluster. 
The red curve shows the X-ray flux limit of the survey, for the median radio 
flux. 
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Fig. 9. — Fraction of radio AGN in cluster cores versus X-ray luminosity. 
The ratio of the number of clusters with radio loud AGN within 250 kpc) to 
the total number of clusters in each X-ray luminosity bin is plotted against 
the X-ray luminosity. To reveal the effect of the redshift-dependent lu- 
minosity cut, the clusters are separated into two subsamples by redshift, 
0.3 < z < 0.5 (blue) and 0.1 < 2 < 0.3. The mean redshift of each 
bin is marked on the plot. Only radio sources with Pi 4 > 3 10^* WHz^^ 
are counted, corresponding to the radio power limit at z ~ 0.5. Uncertainties 
are calculated assuming Poisson statistics. 
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Fig. 10. — Distribution of radio sources around 400SD clusters. The ra- 
dio sources are grouped into two subsamples by redshift of the host cluster: 
0.1 < 2 < 0.3 for the solid hne and 0.3 < 2 < 0.5 for the red dashed line. 
The gray solid line shows the distribution for lower redshift sources, for clus- 
ters with Lx.bol > 10^^ erg s^^ for th e ho st clusters. The numbers have 
been corrected for background sources (J23). Error bars show the standard 
deviations for each point, while those for the gray solid line are omitted for 
clarity. 



spondence between the gray and black curves shows that their 
difference from the radial distribution for the higher redshift 
sample (red dashed line) is not due to the X-ray flux limit. 
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Fig. 1 1. — Distribution of rad io power around 400SD clusters. The sym- 
bols are the same as for Fig.[T()] but here the average power of NVSS radio 
sources within r is plotted against r. Note that the brightest source (3C288, 
with Pi. 4 ~ 700 X 1024WHz-i, seeFig.[3) atz -^ 0.25 is omitted (see 
813.2) . If its radio power were included, the average powers for the lower red- 
shift sources would be increased by 5 X lO^'* WHz^^. The uncertainty in 
the average power is calculated assuming a Gaussian distribution. 

Like Fig. [TOl Fig. [TT] shows the distribution of the average 
radio power around 400SD clusters. The radio power also 
show a modest increase with redshift, but here the impact of 
the X-ray flux limit (see the gray solid line) is more marked. 
Note that the average power for the low redshift sample does 
not include the brightest radio source (3C288, z ~ 0.246), 
which is roughly an order of magnitude brighter than the sec- 
ond brightest source within 250kpc (Fig.Q. If the contribu- 
tion of this source was included, the average power of the low 
redshift sample would exceed that for the high redshift sam- 
ple. We omit this source because the jet power derived from 
its radio power exceeds the jet power determine d fro m its cav- 
itiesby more than an order of magnitude (see il3.3l iLal et al.l 
I2OIOI) . This discrepancy may simply reflect the large intrin- 
sic scatter in the relationship of Eqn. |2] but we are unable 
to pursue this issue here, since our sample lacks a statistically 
significant population of such extremely luminous, short-lived 
sources. 



of isolated ellipticals, which is -- 15%^ Best et"an (l2007h . So 
it is possible that our host galaxies lie at the centers of weak 
cooling flows. This issue will need to be addressed using deep 
X-ray imaging. 

In Fig. [TO] the clusters are divided into two redshift ranges, 
0.1 < z < 0.3 (solid line) and 0.3 < z < 0.5 (red dashed 
line). We find that the average number of radio sources in- 
creases with redshift, consistent with the discussion of Fig.|9] 
in il3.1l Nevertheless, the comparison between the higher red- 
shift and the lower redshift clusters shows that the spatial dis- 
tribution of radio sources around these clusters is not signifi- 
cantly different from that for all lower redshift clusters in the 
400SD sample. In order to show that the redshift dependence 
of the average number of radio sources is not a consequence of 
the X-ray flux limit (i.e., luminosity limit) for the higher red- 
shift clusters, we also plot the distribution of radio sources for 
lower redshift clusters with Lx.boi > 10^'' erg s^^ (gray solid 
line). This X-ray luminosity corresponds approximately to 
the flux limit for the higher redshift sample, so the close corre- 



3.3. Jet power estimated from radio power 

Radio power has been shown to correlate with "cavity 
power," an estimate of the mean power of radio jets based 
on enthalpies and ages of X- ray cavities (iBirzan et al.ll2004 
l2008t ICavagnolo et al.l 120 Iff). Given their close correspon- 
dence to radio lobes, cavities are assumed to be formed 
by AGN outbursts. As a cavity rises buoyantly and ex- 
pands, its enthalpy is expected to be released into the hot 
atmos phere that hosts it (reviewed in iMcNamara & NulsenI 
l2007h . This is one of several mechanisms by which ra- 
dio outbursts may heat the gas. Other mechanisms include 
shock waves or sound waves created by the expansion of cav- 
ities d^rman et al. 2005; M cN amara et al. 2005; F abian et al] 
l2005bnNulsen et al.ll2005allbt iForman et al.l l2007i) and ther- 
mal ization of cosmic rays leaked frorn a radio source 
(Bohr inger & MorfilllfTOSl lMathewsll2009t iGuo & Mathewsl 
i20ili) . Here, we estimate jet powe rs, Pjet, us ing the P1.4 - 
Pjet scaling relation of Cavagnolo e t al.l (120101) . 

log Pjet = 0.75 (±0.14) log P1.4 + 1.91 (±0.18), (2) 
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Fig. 12. — Top: Average jet power (Pjct) of radio galaxies in the 400SD 
clusters. We estimate the jet power of radio s ources from the 1 .4G Hz flux 
in the NVSS catalogue, using the calibration of lCavagnolo et all J20i0). The 
averages are calculated in 5 redshift bins covering 0.1 < z < 0.6, for the 
matched NVSS sources within 250 kpc (black circles) and 500 kpc (red cir- 
cles) of each cluster center For the lowest redshift bin, results are also shown 
for subsamples with Lx > 10^^ ergs^^, plotted as stars. The numbers of 
clusters hosting NVSS sources and the total numbers of clusters in each red- 
shift bin are noted at the bottom/top of the plot. The error bars, estimated 
using a Monte-Carlo method, are dominated by the scatter (cri.4) in the cor- 
relation between jet power and radio luminosity. Bottom: Ratio of jet power 
to X-ray luminosity for 400SD clusters vers us re dshift. The plot shows the 
average jet powers from the top panel of Fig. ll2l divided by the average total 
X-ray luminosity of the clusters in each redshift bin. 
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Fig. 13. — Ratio of jet power to X-ray luminosity versus X-ray luminosity. 
Similai' to Fig |12l but here the sample is grouped by X-ray luminosity. 

which is expected to provide, on average, a minimum esti- 
mate for the total power injected into a cluster by a radio 
AGN. Here, Pjot is in units of lO'''^ erg s^^, and Pi 4 in unit of 
10'*° erg s^ * . The scatter in the correlation between jet power 
and radio luminosi ty is gi.4 = 0.78 d ex. It is dominated by 
the intrinsic scatter jBirzan et al.l2008h . due to a number of ef- 
fects, including aging, and differences in particle content and 
magnetic field. The relationship (Eqn. |2]i is determined over 
6 decades in Pi. 4 (lO^^-lO"*"^ erg s^*), including data for nu- 
clear radio sources for BCGs in cooling core clusters, with Lx 
up to lO**^ ergs^^, and low power sources in giant ellipticals 
with Lx around lO''^ erg s^*. The large scatter in this scaling 
relation implies that the jet power estimate it provides for any 
individual object has an uncertainty approaching a factor of 
10. 

We estimate the jet power for all of the cluster radio sources 
using Eqn. |2] albeit with a large intrinsic uncertainty. Assum- 
ing the errors are random, we expect the uncertainty in the 
ensemble average to be reduced to an acceptable level. More 
importantly, the average jet power calculated for an unbiased 



sample should provide an estimate of the contribution to the 
energy budget of AGN power integrated across the sample. A 
more precise measurement of jet power may be obtained for 
an individual object by analyzing its X-ray cavities. Despite 
being more precise, the measurement is only an instantaneous 
snap shot of its current jet power, and does not necessarily 
reflect the historical average AGN power. The ages of FR 
1 radio AGN outbursts in clusters lie typically between 10^ 
to 10^ y rs, and the lives o f powerful FRII sources are even 
shorter (lO'Dea et al.ll2009i) . These time scales are much less 
than cluster ages, which greatly exceed 10^ years. Therefore, 
an estimate of the time average jet power is required to cal- 
culate the total AGN energy input to the clusters over cosmic 
time. This quantity may be estimated by averaging jet power 
over a sufficiently large and unbiased sample of clusters, as 
we have done here. 

In Fig. [T2I we show our estimate of the average jet power 
(Pjot), using a Monte-Carlo method that allows for the uncer- 
tainty in the radio flux and the parameters in Eqn.|2l the distri- 
bution of radio spectral indicefl and the large intrinsic scatter 
(o'i.4) in the relation (Eqn. |2]i to estimate confidence ranges. 
Overall, the confidence ranges shown in Fig.[T2]are dominated 
by CT1.4. Note that we use the arith metic mean value of Pjct, 
excluding the brightest source (see i l3.2b . While this gives a 
mean power that is substantially greater than the "mean" of 
the log-normal distribution of jet powers, assumed implicitly 
in the correlation (Eqn.|2]i, our purpose is to determine the av- 
erage energy deposited by AGN, not its distribution, so that 
the arithmetic mean is the appropriate estimator In order to 
quantify the potential underestimate, we calculated the mean 
jet power excluding the 10% most powerful radio sources, and 
found that the number drops by 5%. 

We find that the average jet power in the cores of the 400SD 
clusters is 2x 10'*"'ergs^*anddoesnot showasignificantred- 
shift evolution. The slight increase in the jet power at higher 
redshift reflects the increased AGN fraction shown in Fig. |9] 
and Fig. [TO] The typical total X-ray luminosity of the 400SD 
clusters (Fig. [U 0.2 < z < 0.6) is - 3 x 10^"* ergs-\ so 
that the average jet power in these clusters is about 50% of 
their X-ray luminosity (Fig.[T2]i. Comparing the two plots in 
Fig. [T2I we find that the higher ratio of jet power to X-ray 
luminosity for the lowest redshift bin shown in the bottom 
panel of Fig.[T2]is caused primarily by selection effects. Us- 
ing an X-ray flux limit reduces the average X-ray luminosity 
at lower redshifts, while lower luminosity clusters at higher 
redshifts fall below the detection threshold. Excluding low 
luminosity clusters (Lx < 10^"* erg s^*) for the low redshift 
bin 0.1 < z < 0.2 gives an average jet power (star symbols 
in Fig.[T2ll that is comparable to the values at higher redshifts, 
consistent with the average radio pow er be ing largely inde- 
pendent of the redshift, as argued in i i3.2l (see Fig. [TTI ). It 
suggests that the average jet power does not depend on the 
X-ray luminosity of a cluster If so, the impact of AGN on 
the gas will be greater in lower luminosity, hence less mas- 
sive, systems. Fig.fljjshows that the minimum jet power may 
exceed the power radiated by the ICM in group-scale objects 
withLx < 10*** ergs^*. If the jet power is mostly retained by 
the hot atmospheres and not radiated away, AGN can cause a 
net increase of entropy in low mass systems. 



3.4. AGN heating 



^ A Gaussian distribution with a mean of 0.75 and a it of 0.10 <Condonl 
119921) is assumed for the spectral indices. 
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TABLE 1 
Energy per particle 



Nci 


-f'X.bol 
lO'^'ergs-i 


Ndct 


<250kpc 
heating 
keV par^ ^ 


Ndct 


<500kpc 
heating 
keV par^ ^ 


65 

57 
47 


0.47 
1.73 
6.97 


19 
19 
16 


0.23 ±0.12 
0.12 ±0.06 
0.17 ±0.10 


56 

42 
36 


0.39 ±0.17 
0.27 ±0.13 
0.27 ±0.14 


169 


2.71 


54 


0.17 ±0.08 


134 


0.30 ±0.14 



From the jet power, we can estimate the mean total energy 
per particle injected by the radio sources (Table [T]i 



Ejct = .Pjct iz/.<"lp/Afga 



(3) 



where n = 0.59 is the mean molecular weight and Wp is the 
proton mass. Pjct and iiigas are the mean jet power and gas 
mass, respectively, for the cluster^ The gas mass is calcu- 
lated fro m the X-ray luminosity assuming the LX-M500 re- 
lation in IVikhlinin et al.l (120091) and a gas mass fraction of 
0.12. We excluded clusters beyond z = 0.6 because their 
numbers are small. The duration of energy injection, t^, is 
estimated as the time interval from z = 0.6 to 0.1 (4.4 bil- 
lion years). As might be expected from the discussion of 
Fig. [T3] Ejct decreases with X-ray luminosity, although the 
trend is weaker due to the luminosity dependence of the Lx- 
M500 relation. AGN heating from within the core alone can 
reach about 0.2keV/particle for poor clusters with typical X- 
ray luminosities of 5 x lO^'^ergs^^. The number increases 
to ^ 0.4keV/particle if the aperture increases to 500 kpc. We 
have also combined the data for the three X-ray luminosity 
bins, giving the average AGN energy input per particle for the 
entire sample in the last column of the table. 

We expect the average energy injected by radio sources to 
be underestimated somewhat due to the radio luminosity cut, 
particularly for the high redshift clusters. We examine this 
effect by calculating the contribution from the weaker radio 
sources (Pi. 4 < lO'^'* W/Hz) in clusters lying between red- 
shift 0.1 and 0.3. We further assume that the radio luminosity 
function does not evolve with redshift. Including the contri- 
bution of the weaker sources, the average energy injected by 
radio sources into clusters at 0.3 < z < 0.6 is boosted by 
7% within 250 kpc radius and by 11% within 500 kpc. This 
correction ignores radio sources lying below our 3 mJy radio 
flux limit. Lowering the limit to 2 mJy and weighting to off- 
set to the flux-dependent completeness of the NVSS survey 
([Condon et al. 1998), the increase in the average injected en- 
ergy is negligible (< 2%). Since the dependence of jet power 
on radio power according t o Cavagnolo' s calibr ation is steeper 
than found previously by iBirzan et al.l (120081) . the contribu- 
tion of weak source is expected to be less important than was 
found, for example, bv lHart et al.l (|2009). 

We note that the numbers in the table are calculated us- 
ing the mass of gas within Rsoo- The radio sources consid- 
ered here deposit most of their energy directly into the clus- 
ter cores, but current ignorance of the significant heating pro- 
cesses means that we do not know where the bulk of the en- 
ergy will ultimately reside. If it did remain within the cluster 

' The mean jet power we estimated is defined within an aperture projected 
on the sky, but the gas m ass here is defined in three dimension. From Fig. O 
and discussions in q3.2l the radio sources in clusters are dominated by central 
sources. Therefore, the difference between projected and three dimensional 
mean jet power should not be significant, especially for small aperture like 
250 kpc. 



cores, the energy deposited per unit mass would be approxi- 
mately ^ 5 time^3) larger than the values given in the table, 
with a commensurate increase in its local impact. However, 
our primary objective is to demonstrate that the energy input is 
significant, even if the energy ends up being spread uniformly 
throughout much of the cluster atmospheres. 

We have limited our jet power calculation in the redshift 
interval of the sample (0.1 < z < 0.6). In principle, the 
energy injected by radio sources shoul d be traced back to the 
time w hen BCGs formed (z r-^ 2.0, e.g. Ivan Dokkum & Franxl 
1200 II) . Ignoring AGN evolution, the energy injected per par- 
ticle since z = 2 would be at least twice the values listed in 
Table [T] This is a conservative estimate, si nce the fraction 
of powerful AGN increas es with redshift (e.g. lGalametz et all 
l2009HMartini et al.ll2009h and the contribution of AGN heat- 
ing is expected to be significantly greater in the past. Further- 
more, we should take into account the evolution of clusters 
over this time scale, so that we cannot simply scale the val- 
ues in the first three rows of Table [T] which are derived for 
fixed X-ray luminosity ranges. Naively, if we assume that 
the mass function of 400SD clusters is representative of the 
"integrated" mass function of clusters since z = 2, we can 
extrapolate the numbers in the last row of Table [1] to provide 
a rough estimate of the minimum AGN energy injected per 
particle since z — 2. The values we find are 0.4keV/particle 
fr om AGN \ yifliin 250 kpc. 

IWu et al.l (12000 ) found that the minimum excess energy re- 
quired to break self-similarity is ~ IkeV/particle. There- 
fore, our results suggest that the continuous input of energy 
from AGN activity at the current rate over the ages of clusters 
would provide a significant fraction (40%) of the excess en- 
tropy (heat) found in the hot atmospheres of clusters. It also 
implies that AGN heating may be a factor leading to the dearth 
of strong cooling flows in distant clusters (e.g. IVikhlinin et all 
2007; Samu ele et al. 2011) . It is unclear whether this heating 
is being supplied by a self-regulated feedback loop. 

4. SUMMARY 

We have correlated the positions of clusters from the 400SD 
cluster survey lying between redshifts 0.1 and 0.6 against ra- 
dio sources with fluxes above 3mJy from the NVSS. Radio 
sources are detected within 250 kpc of the center for 30% of 
the clusters and within IMpc for 50 to 80%, depending on red- 
shift. The first valu e agrees wi t h the fraction of radio sources 
in BCGs found by iBest et al.1 (I2007h for SDSS clusters, but 
lies well below the nearly 1 00% detection rate in stron g cool- 
ing flows (e.g. Burns! 1 1990; Edwards et al. 2007; Dunn et ajj 
2006; Mittal et al. 2009; Cavagnolo et al. 2009). We found 
that radio power of radio sources in the 400SD clusters does 
not correlat e with the cluster X-ra y lum inosity. For compar- 
ison, Mittal et all (|2009) and Sun (^2009") found strong corre- 
lations between AGN radio power of BCG host galaxies and 
the X-ray luminosities of their bright cooling cores, and the 
correlation disappeared in systems with weak or absent cool- 
ing cores. Both of these facts support our conclusion that the 
400SD sample is composed primarily of weak or non-cooling 
core clusters. 

We have found no significant correlation between the in- 
cidence of radio sources in clusters and the X-ray luminosity. 
Since the X-ray luminosity of a cluster correlates with its rich- 
ness, this implies that richer clusters have fewer radio AGN 

'" The factor is estimated from the ratio of mass within R500 and 250 kpc 
radius assuming a simple density profile p oc r~^. 
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per galaxy. In addition, the density profile of radio sources 
in clusters peak at the cluster center, implying that BCGs are 
the dominant radio galaxies in clusters. Apart from the BCG, 
the rich cluster environment is hostile to the formation of ra- 
dio AGN, so that the number of radio AGNs does not increase 
with the number of galaxies. 

We have found only modest correlation between radio 
power and incidence of radio detection in 400SD clusters with 
redshift. Although it is difficult to completely rule out bias 
in a flux-limited sample, our results are qualitatively consis- 
tent with redshift evo lution of the radi o AGN detection rate in 
clusters by e.g. Gala metz et alJ (120091) . 

We have estimated the average jet power and the rate of 
AGN heating in clusters using the scaling relation of Eqn. |2l 
The average jet power is about 2 x 10*"'ergs^^for radio 
AGNs within 250 kpc radius of the 400SD clusters and about 
4 X IC^^ergs^^for sources within 500 kpc radius. We found 
that the average jet power and AGN heating rate do not corre- 
late with total X-ray luminosity of the 400SD clusters. There- 
fore, the heating rate per particle will be larger in less massive 
systems. The AGN heating within the core of clusters can 
reach 0.2keV/particle for poor cluster with typical X-ray lu- 
minosities of 5 X lO''^'^ erg s^^. These numbers are calculated 
within the redshift range of our sample, i.e. 0.1 < z < 0.6. If 



we extrapolate the result to redshift z = 2 ignoring AGN evo- 
lution, the integrated AGN heating per particle will increase 
by a factor of 2. 

If the heated gas is unable to cool quickly enough, the en- 
tropy of their hot atmospheres will rise above the values ex- 
pected from gravitational heating alone. In fact, we found in 
i l3.4l that the amount of AGN heating of the hot atmospheres 
is a significant fraction (40%) of the he ating requi red to "pre- 
heat" clusters (lKaisedll991l; IWu et alJ lTOOO; McC arthy et al] 
l2002h . This so-called preheating phase is thought to occur 
durin g the ep och of galaxy formation at redshifts of 3 and be- 
yond (Kais3[i99il)- Our results suggest that the AGN heat- 
ing of cluster atmospheres long after the epochs of galaxy 
and cluster formation, and throughout the formation history 
of clusters are significant. Thus, the heating that apparently 
broke the self-similarity of cluster scaling relations appears 
to have occurred continuously, and not necessarily at a single 
epoch. 
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